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^^ Abstract 

en ' 

0^ \ Polarized beam ^p Compton scattering provides a theoretically clean way to 

Q \ extract the isovector parity violating pion-nucleon coupling constant h^ji^^- 

^D ' This channel is more tractable experimentally than the recently proposed ex- 

traction of h\ljj^ from the Bedaque-Savage process — polarized target 773* 
compton scattering. The leading parity violating effect is calculated using 
O \ Heavy Baryon Chiral Perturbation Theory {HBxPT). The size of the asym- 

^ I metry is estimated to be ~ 4 x 10^^ for 120 MeV photon energy. 

> 
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The isovector parity violating (PV) pion-nucleon coupling constant /i^^^y is responsible 
for the longest range part of the A/ = 1 P V NN forces [|T]-^ • It is expected to give dominant 
contributions to low energy quantities such as nucleonQ and nuclear anapole moment |]6|-[ll 



(1) 



and PV neutron radiative capture np — >■ dj. However, past attempts to extract hl^lf^ ^^^ ^^^ 
satisfactory (see 0JT^,[TB| reviews). In many-body systems, several PV effects are enhanced 
and have been detected. On the other hand, the theoretical analysis is complicated. The 
extractions from ^^F [[T4| , p!5| and ^^^Cs |[16HT9| systems differ by an order of magnitude with 
large uncertainties, while the measurement in the ^°^T/ system gives a null result |20[. In 



fewer body systems, the theory is more under control but the PV effect is smaller such that 
previous measurements could not reach the required precision 



mm 



However there are 

several high precision new measurements under preparation or execution including ftp -^ d'j 
at LANSCE [^], 'yd -^ np at JLab [^, and the rotation of polarized neutrons in helium at 
NIST. It is expected that these experiments will put tight constraints on the value of /i^ 



NN- 

In the single nucleon system, a new PV observable was recently suggested by Bedaque 
and Savage ^7j. They found that the polarized target •y^ —>■ 7p Compton scattering 
asymmetry measurement (calculated to be ~ 5 x 10"*^ for 100 MeV photon energy) would 
determine h^j^j^ with an estimated 15% uncertainty. To control systematic errors, the 
difference in cross section for the proton spin polarized parallel and antiparallel to the 
direction of the incident photon must be measured during a short period of time with 
only the target polarization direction changed. To achieve this, the rapid flipping of the 
target polarization should be employed. This is unpractical with the currently available 
experimental techniques and polarized proton targets, thus the Bedaque-Savage (BS) process 
is not favored experimentally. On the other hand, rapid flipping of beam helicity is a standard 
technique already employed in many parity violating experiments. Thus the polarized beam 
experiment, ^p -^ 7p, is experimentally more tractable [^ 



.(1) 



Given the great interest in the 
determination of h)^j^^^ we investigate this ^p parity violating process using Heavy Baryon 



Chiral Perturbation Theory (HBxPT) il,|33|. 

We start with reviewing the symmetry constraints on the Compton scattering process 



7(A;,e) +p -^ 7 (A;', e') +p 



(1) 



where [k, e) and [k', e') are the initial and final photon momenta and polarization vectors. It 
has been known for a long time that there are ten time reversal invariant structure functions 
in the transition amplitude for this process. Six of them conserve parity while four of 
them violate parity. This can be seen easily by the following exercise on helicity amplitude 
counting. Using |Ai,A2) to denote a state with photon and proton helicity Ai and A2 in 
the center of mass frame, sixteen helicity amplitudes (A'^, A2IA1, A2) can be constructed for 
proton Compton scattering. These amplitudes transform under time reversal as 



(Ai, A2IA1, A2) — > (Ai, A2IA1, A2) 



(2) 



^ The isoscalar nucleon anapole moment is dominated by h\^lj^, but not the isovector anapole mo- 



,(1) 



ment which is more directly relevant for the SAMPLE electron-proton 
PV experiments. 



and electron-deuteron 



which is equivalent to taking a transpose transformation of the 4x4 matrix. Thus ten time 
reversal invariant amplitudes can be constructed through the linear combination 



Mt{\i, \2, K, K) = (A'l, A'^lAi, A2) + (Ai, A2IA;, A'2) 



(3) 



One can further separate these into PC and PV amplitudes. Two of the ampli- 
tudes, Mt{+1, +1/2, -1, -1/2) and Mt{+^, -1/2, -1, +1/2), have an additional symme- 
try. They are invariant under parity transformation. 



(Ai, A2IA1, A2) -^ {—^i, — A2I — Ai, — A2) 



(4) 



after having been made to conserve time reversal invariance. The other eight amplitudes 
can be grouped into four PC and four PV amplitudes using similar linear combinations to 
that of eq.(^. This demonstrates the well-known result that there are six independent PC 
and four independent PV structure functions satisfying time-reversal invariance in a proton 
Compton scattering process. 

In the center of mass frame, the six PC structure functions can be chosen as 



rppc ^ JY 



Ai e-e'* + A2k • e'* k' • e + iAacr ■ (e'* x e) + iAicr • ( k' x k)e ■ e' 



+i^5cr ■ 
+iA6a 



y* X k)e ■ k' - (e X k')e'* ■ k 
(e'* X k')e ■ k' - (e X k)e'* ■ k 



A^ 



(5) 



where N is the proton spinor, a is the Pauli matrix acting on the nucleon spin index, k and 
k' are the unit vectors in the k and k' directions, and the Coulomb gauge (eo = eg = 0) is 
used. The PV structure functions can be chosen as 



-J^sk ■ e'* k' ■ e a • (k + k') - ij^^e x e'* ■ (k + k') 



N 



(6) 



The JFi_3 structures were first given in ref. ||2^. The interference between Ai^2 and J-'i.s 
contributes to the BS process. For the polarized beam process ^p —>■ ■jp considered here, 
the contributions are from the interference between Ai^2 and JF4 and between ^3.6 and J-'i^. 

Since we are interested in the low energy behavior of the proton Compton scattering 
process, chiral perturbation theory provides a natural framework with which to work. As a 
low energy effective field theory of QCD, chiral perturbation theory captures the symmetries 
of QCD and describes low energy observables by derivative and Chiral expansions. The 
SU{2)l X f/(l) symmetry structure of electroweak interactions can also be incorporated 
with the weak boson exchange described by contact interactions while keeping the photon 
as dynamical degrees of freedom in the Chiral Lagrangian. 

HBxPT has been applied to the calculation of several Compton scattering observables. 
The PC structure functions have been calculated up to next-to-leading order (NLO), 0{e^p) 
PB| , and are listed in the Appendix. The proton Thompson term is the only contribution at 
leading order (LO) and contributes to Ai only. The NLO contributions come from the tree 
level diagrams, pion-nucleon loop diagrams and the Wess-Zumino term. They contribute 








FIG. 1. The leading order contribution to parity violating sturcture functions T\-J-j, in ^p 
Compton scattering. The solid square is the weak operator with coefficient h^j^j^- Wavy lines are 
photons, solid lines are nucleons, and dashed lines are pions. The crossed graphs are not shown. 
Graphs with photons from the strong vertex, or insertion of the two-photon-pion vertex vanish in 
the v.A = gauge, and thus are not shown here. 
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FIG. 2. The first non-vanishing order (NLO) contribution to parity violating sturcture functions 
T^ in 7p Compton scattering. The features of the graphs are as defined in fig^. The photon nucleon 
couplings are magntic couplings. 

to all the PC structure functions ^i-e- These structure functions can be determined exper- 
imentally (see the Appendix for more details) thus the uncertainty from the PC part can 
be eliminated completely. Here we use the HBxPT result only for the sake of estimation 
of the asymmetry. For the PV structure functions, the LO {OlGpe"^) with Gp the Fermi 
coupling constant) contributions have been calculated in ref. |2^. We also list them in the 
Appendix. They arise from the pion loop diagrams shown in fig.0 and contribute to Tx-T-^. 
^4 is an additional quantity we will need to compute. Its contribution comes from the pion 
loop diagrams shown in fig.|^ and the effect starts at NLO. Now we give some details of 
computing the T^^'~' . 

The PC part of the relevant Lagrangian is 

1 2 

+ ^AT {{v ■ Df -D']N- -^N [S^, 51 [/.o + f^ir,] NF,^ + ■ ■ ■ (7) 

where the pion decay constant F^ = 93 MeV, pion-nucleon coupling constant qa = 1.26, 
A^ is the isospin doublet of the nucleon fields with velocity v, D is the covariant derivative 
with gauge coupling on the proton as D^^N = [d — ieA)^N, and S is the covariant nucleon 
polarization vector. In the proton rest frame, v^ = (1,0,0,0), S^ = (0,cr/2). /io = (/ip + 



/i„)/2 and fii = {fip — /U„)/2 are the isoscalar and isovector magnetic moments in nuclear 
magnetons, with fip = 2.79 and /i„ = —1.91. The elhpses denote terms with more pion fields 
and insertions of higher powers of derivative and pion mass. Massive hadronic excitations 
such as kaons and deltas are "integrated out". Their effects are encoded in the higher 
dimensional operators. 

The non-leptonic PV part of the relevant Lagrangian is 

j^pv ^ 'biEKe^^^T.^r^N + ■■■ = -th'-^l^n+p^n + h.c. + ■■■ (8) 

v2 

where the ellipses denote terms with more pion fields and derivatives. This Lagrangian was 
first given in ref. [Q with a different phase convention for the pion field. We adopt the same 
convention as refs. [|l^,|Tl|]. /i^/rTv ^^^ estimated by matching onto four quark Fermi theory 



and was found to be dominated by s quark contributions, 



"'wNN 



^^ ~ 5 X 10- , (9) 



where A^ ~ 1 GeV is the chiral perturbation scale. This estimation is consistent with the 



"best value" obtained in ref. M and close to one result PS| from QCD sum rules. A recent 



calculation in the SU(3) Skyrme model yields /i^^jv^ ~ 0.8-1.3 xlO ^ [^. The radiative 



W „ , n s 1 Q NX 1 n-7 



,(1) 



correction on the /ij^^v^ vertex is discussed in ref. [00 



The /i^^Tv term is the only term in the C^^ with a non-derivative pion-nucleon coupling; 
it is expected to dominate over the other contributions to JF4. The pion loop diagrams in 
fig.| gives 



^4 



NLO 



8V2mNF^ 



TC^ 



ml I . _^ uj 



uj sm 

uj \ m. 



(10) 



with a magnetic photon-nucleon coupling, uj is the photon energy in the center-of-mass 
frame. 

In this calculation delta contributions are encoded in the higher-order operators and will 

contribute through higher-order diagrams. If one considers the delta-nucleon mass difference 

A ~ 300 MeV as a light scale as itIt^ and uo are, then one needs to sum factors of itlt^/ A and 

uo/ A to all orders. This can be done by including delta as a dynamical degree of freedom 

32,34]. In this expansion, the delta diagrams will contribute to J-'f^Jf^ through vrA loop 



diagrams and tree diagrams with unknown 7A^A couplings. In the expansion with which we 
work, the A is considered as a large or heavy scale, so factors of m.^j IS. and cj/A are treated 
perturbatively. Thus below pion production threshold [oj < m.„), the delta would contribute 
a factor of (r7i^/A^,co'^/A^) ~ 25% correction to J-'^^'^. This is the dominant source of the 
uncertainty. 

The PV asymmetry can be defined by the difference in the cross section (in the center- 
of-mass frame) for photon helicity A^ = -|-1 and —1 normalized to the sum 



da ,^ , da 






dH^'^' = +1) + dn*"- 



Again, the PC structure functions Ai.q can be extracted from experiments to further reduce 
theoretical input and ehminate the uncertainties from the PC part. Here we plug in the 
PC HBxPT result in order to get an estimation of the size of the asymmetry. In this 
treatment, the helicity asymmetry contribution starts at NLO, 



cos 



A,, {uj, 9) = 1^^ Re {^fo^f ^o* + <^° [^i^«* - 2 J-2^°* - J-3^«* ( 

\A^o\\l + cos^ 9) ^ ' ^ 3 L 1 3 V 

+^™ [^1^°* (1 + cos^) + J-2^°' (1 - 3 cos^) + J^i""' (1 - cos^ 9)] 

-A^"-^ f^f'^* (1 + COS0) + J^^^' (3 - cos^) + J^^^' (1 - cos^ 9)] 

-<^'^^2^^*(l + COS^)} 

4(1 + COS^) ^^ ^ J^^^LOjrLO'' ^ ^^0_^7VL0*| _ (^2) 



i^for(i+cos2 

We have used the LO result for the PC cross section, so 



da ,^ , da 



'-1 



|2 ,. 9^, O? 



^j^KI (l + cos^«) = ^(l + cos-^ 



(13) 



It is instructive to study the low energy limit [lo <^ vn.^^ of the asymmetry. Keeping the 
first term in the uo/m.,^ expansion, 

^f^ = -^ , ^f-^~[l + 2K,-(l + M^cos^]^ , 
anlo _ _ anlo (1 + Up) e LO anlo (1 + Hp) e uj 

2 uW 2 2 uW 3 

•^1 = -^9 ~ ;= 1 '^d ~ H ^ 1 (-'-4) 

2AV2n^F^ml ^ 24V27r2F^Mjvm2 ^ ^ 

with J-'^''^* = 0{u}^ /m^) and Kp = fip — 1. The inverse power of niT, dependence in the jFs 
explains why there are no intrinsic unknown two-photon-two-nucleon counterterms at this 
order. In this low energy limit, the asymmetry has a simple form 

24v 27r^i'^m^(l + cos"' 9) 

The vanishing of A^^ at backward angle {9 = vr) is a consequence of time-reversal invariance 
and is a general property for all values of the photon energy. For back scattering, the change 
of photon spin direction corresponds to a A J = 2 operation and hence, is a forbidden tran- 
sition for the proton matrix element. Thus both photon and proton spin directions do not 
change but the helicities change signs. Using eq.(0) for the time-reversal invariant amplitude 
and (A']^,A2) = (— Ai,— A2), this amplitude conserves parity and does not contribute to the 
PV asymmetry. 
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FIG. 3. The estimated photon helicity asymmetry A^^ defined in eq.(ll) shown as a function 
of photon energy with the photon reflection angle 9 = tx 12 in the center-of-mass-frame. The naively 

(1) —7 

estivaated value of h^j^j^ = 5 x 10 is taken as input, and a HBxPT estimation used for the parity 
conserving amplitude. 



For a numerical estimation of the magnitude of the asymmetry, we consider 6 
where 



A 



"'ttNN 



(u « m., -) ~ -8.8 X 10-9 JhiNIL^ 
\ '2/ I 5x10-7 



to \3 

KTOMev) 



n/2, 



(17) 



with ~ 25% uncertainty. 

In fig.^, we show the photon energy dependence of the estimated asymmetry at 
6 = 7r/2. Assuming the naive size for hlj^j^ estimated in eq.(|), the asymmetry is 

A^^ (120 MeV,7r/2) 3.8 x 10-^ with the higher-order uncertainty ^(m^/A^wVA^) ~ 

25%. Note that very near the pion production threshold, resummation of terms with powers 
of 



mz 



2mN iuo 



m^ 



(18) 



is required in order to shift the pion production threshold from r/i^ to m^r + 



m: 



2m 



fin the 



N 



laboratory frame) to recover the recoil effect. The resummation procedure is well known [|35| . 
Without this resummation, we should restrict ourselves to cj < 130 MeV such that the factor 
in ([T8|) is sufficiently less than 1. To probe how far away from the threshold our calculation is 
still under control, we compare the results of expanding the A^^ to NLO to that of expanding 
the total amphtude to NLO, then square it to get the A^^ (the A2 J^4 interference term is 
included in the latter case). If the difference between these two quantities is consistent with 
the estimated higher order contribution, then the result without extra resummation will be 
valid. We find the difference increases with u and reaches 14% at cu = 120 MeV, so our 
expansion is presumably still useful up to 120 MeV. 
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FIG. 4. The angular distribution of the estimated photon helicity asymmetry A^^ defined in 
eq.^T\) and calculated in HBxPT in the center- of -mass- frame for 100 MeV photon energy. The 

(1) —7 

naively estimated value of h^^j^ = 5 x 10 is taken as input. 



In fig.^ we show the angular distribution at cu = 120 MeV for the same estimated value 
of h^j^j^- The maximum asymmetry is near 9 = n/2 but slightly biased towards the forward 
direction. The asymmetry vanishes at = tt as required by time reversal invariance. 

In conclusion, parity violating ^p -^ 7p Compton scattering provides a theoretically 
clean way to extract h^j^j^. The dominating source of the PV effect comes from the PV pion 
loop contributions. The magnitude of the helicity asymmetry is estimated to be ~ 4 x 10~^ at 
120 MeV photon energy with ~25% uncertainty for a natural size h^^lfj^, under the framework 
of HBxPT. Thus we have found a model independent way to constrain h\^j^j^ with ~25% 
uncertainty. We note that the ~25% uncertainty is dominantly due to the delta and can in 
principle be reduced by inclusion of the delta as an explicit degree of freedom. Unfortunately 
this would require additional experiments to measure the PV 7A^A coupling. However, even 
with ~25% uncertainties, ^p Compton scattering will greatly improve our understanding of 
^ttNN- ^^ ^^6 optimistic that this experiment is feasible for current experimental techniques 
and facilities. 
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Appendix 

The following PC structure functions computed to 0{e^p) are taken from eqs. (4.28a-g) 
of ref. IQ (with a typo in eq. (4.28b) corrected). 



A{ 



LO 



M, 



N 



jNLO _ 9Ae 



^- A/m2 -UJ^ + 



2Ml-t 



1 V-t 
- arctan 

2 2m^ 



l-zW-t 



— / (iz arctan , 

/o 2v/M^ 



LU'^Z'^ 



(19) 



^: 



7VLO _ ^ , e^^i^^ t - 2m^ 



(iz 



arctan 



(l-2)v^ 2(l-;z)v/t(ZJ2^ 



m:f 



2v/m2-a;2z2 4m2 - 4^ V - t(l - 2^1^^ 



A 



NLO 



[l + 2/?^ 



'1 + Kp) COS 



2M| 

+ ^!^^4gij^2^ [ dx [ dz 



e^gAtuj 



+ 



e^9A 



47r2F2 







87r2F2(m2-t) 87r2F2 
x(l — x)z(l — z)^ 



Ml . 2 ^ 
arcsm uo 







W^ 



. ujz ijjzW 



(21) 



^; 



TVLO 



e2(l + K„)2^ 



2M2, 



+ 



47r2F2 



ds / dz- 



z{\ — z) UJZ 



(22) 



A 



NLO 






e^^A^^ 



+ 



e^gW 



87r2F2(m2-t) 8772^2 



1 .1 

da; / dz 
Jo 



„ 9 ^xil — x)z{l — z)^ f . UJZ ujzW 

+2^ cos^^ ^^^ ^ (^arcsm ^ + -^ 



'l-z) 



ijJZ 



W 



■ arcsm 



R 



(23) 



A 



NLO 



e^uj 



2M% 



:i + ^^o) + 



e^QAOj'^ 



^'M£f\. 



-2uj' 



8TT^F^{ml - t) 87r2F2 j^ 

xil — x)z{l — z)'^ ( . ijJZ ujzW 



dz 



W^ 






^aresm- + — 



— Z) . lOZ 



(24) 



with 



j\i 



t= {k- k'Y = -2iu^ (1 - cos^) , 
W = ^Jml - uj'^z'^ + t(l - zfx{x - 1), R = ^/ml + t{l-zyx{x-l). 



(25) 



Expressions of PC Compton scattering for the unpolarized differential cross section, 
da/dQ, proton-photon spin parallel asymmetry A\\, and proton-photon spin perpendicular 
asymmetry A± are given in terms of PC structure functions in eqs. (4. 18) and (4.19) of ref. 
These structure functions can be extracted experimentally. For example, below pion 



production threshold {u < m.„), one can extract Ai and ^3 by measuring da/dVt and A\\ 
at ^ = 0. Since Ai are all real below pion production threshold and A\\ oc ^^ + ^3 and 
A±_ oc AiA^, measurements of ^11 and A±_ are sufficient to extract Ai and ^3. 






The following PV structure functions are given by ref. |27 



^ , dx 

4v^7r2F, Jo JO 

„2„ 7,(1) rl 

e yAii-nNN 



e gAh^NN 



dx 



dx 



l-x 



dy (1 — 2y)uj [X(— 1; xu, w?) — X(— 1; — xtu, rn?)\ 
dy y uj [X(— 1; xuj^ m^) — X(— 1; —xcu, fh? 



l-x 



dy y {1 — X — y) {2y — l)u;^ [X(— 2; xuj, rn^) — X(— 2; —xuj, m^ 



2v^7r2F, JO JO 

rr? = ml + 2y(l — x — y) to'^ (1 — cos 9) , 

where the functions X(a; b, c) are defined by Jenkins and Manohar in ref. ||32 | 

I{a;b,c) = dX (A2 + 2A6 + c)" 

Jo 

X(-l;A,m2 



X(-2;A,m2) 



1 /A - \/A2 -m^ + ie 

2VA2 -m-^ + K ^^\A + VA'^-m'^ + ie 

A 



2 (A^ — m? + ie) ym,^ 



— -X(-l;A,m^ 



(26) 



(27) 
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